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ON  A  NONLINEAR  THEORY  OF  ELASTIC  SHELLS 


By 


W.  L.  Halnwrlght 


Summary 

This  paper  is  concerned  with  a  nonlinear  theory  of  elastic  shells 
with  small  deformations  whose  material  response  is  nonlinear.  The  develop' 
ments  are  carried  out  under  the  Love-Kirchhoff  hypothesis.  General  consti¬ 
tutive  equations  are  derived  in  which  the  geometrical  properties  (due  to 
defonsation)  and  the  material  characteristics  are  separable.  Through  this 
separability,  it  is  shown  how  to  extract  constitutive  equations  of  pre¬ 
determined  types.  Particular  examples  are  followed  by  a  discussion  of  the 
membrane  theory. 


1.  INTRODUCTION 

In  a  recent  paper,  Zema  [i960]  dealt  with  a  nonlinear  theory  of 
elastic  shells  which  has  as  its  basis  the  proposal  of  conatitutive  equa¬ 
tions  which  are  at  once  geometrically  linear  and  physically  nonlinear.  The 
work  of  Zerna  [i960],  however,  is  somewhat  limited  in  scope  since,  in 
particular,  it  contains  only  a  third  order  correction  to  the  linear  theory 
with  no  indication  of  how  a  second  order  or  higher  order  corrections  are 
to  be  obtained.  In  addition,  Zerna's  constitutive  equations  upon  complete 
linearization  reduce  to  a  version  of  those  commonly  referred  to  as  Love's 
first  approximation  which,  as  brought  out  in  more  recent  investigations  of 
Naghdi  [1962],  do  not  satisfy  all  requisite  invariance  requirements. 
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In  Che  present  paper,  through  the  artifice  of  a  series  expansion  for 
Che  strain  energy  function,  we  derive  in  Section  3  general  nonlinear  con¬ 
stitutive  equations  for  isotropic  elastic  shells  whose  mechanical  behavior  is 
(in  the  sense  of  Zerna)  geometrically  linear  and  physically  nonlinear.  The 
expansion  of  the  strain  energy  function  just  mentioned  is  effected  with  the 
aid  of  certain  results  deduced  in  the  Appendix  (Section  6),  where  an 
expression  for  a  potential  is  derived  relating  two  3^3  symmetric  matrices. 
The  notation  employed  throughout  the  paper  is  essentially  that  of  the  paper 
by  Naghdi  [I962]  mentioned  above. 

Specifically,  Che  content  of  the  paper  is  as  follows:  Certain  pre- 
limlnary  results  which  are  needed  subsequently  are  given  in  Section  2. 

As  mentioned  above,  Section  3  contains  the  derivation  of  the  general  con¬ 
stitutive  equations.  Here,  in  addition  to  invoking  the  Love-Kirchhof f 
hypothesis,  it  is  assumed  that  Che  displacements  and  their  gradients  are 
small  but  compatible  with  this  the  constitutive  equations  are  fully  non¬ 
linear.  Also  included  in  Section  3  is  an  investigation  of  the  invariance 
requirements  of  Che  general  constitutive  equations.  In  Section  U,  we  first 
discuss  a  systematic  procedure  for  determining  approximate  constitutive 
equations  of  various  order  in  a  manner  which  results  in  separation  (in 
product  form)  of  the  purely  geometrical  effects  (due  to  deformation)  and 
the  material  response  of  the  shell.  Comparison  is  made  with  the  results 
of  Zerna  [I960]  and  Che  special  cases  of  Che  first  and  second  order  theories 
are  discussed;  the  first  order  (or  linear)  theory  yields  the  Flugge-Lur 'e- 
Byrne  equations.  In  Section  3,  we  specialise  the  nonlinear  constitutive 
equations  to  those  appropriate  for  the  membrance  theory  of  shells.  He  further 
show  that  for  static  problems  of  the  membrane  theory,  even  when  the  con¬ 
stitutive  equations  are  (physically)  nonlinear,  the  determination  of  dis¬ 
placements  may  be  reduced  to  the  solution  of  a  system  of  linear  differential 
equation. 

In  deriving  the  general  constitutive  equations,  the  use  of  power  series 
form  of  the  strain  energy  function  offers  certain  advantages  which  should 
be  elaborated  upon.  Firstly,  in  order  to  obtain  the  ( two-diisens tonal) 
shell  constitutive  equations  from  the  appropriate  strain  energy  function  in 
3-space,  the  power  series  representation  permits  an  explicit  integration 
across  the  thickness  of  the  shell.  Secondly,  through  this  procedure,  there 
results  a  separation  between  the  geometrical  properties  (due  to  deformation) 
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and  those  representing  the  material  response  so  that  desired  geometrical 
and  physical  approximations  may  be  introduced  independently  of  each  other. 
Also,  it  should  be  noted  that  the  power  series  representation  of  the  strain 
energy  function  (discussed  in  general  in  Section  6)  is  independent  of  the 
present  application  and  may  be  used  in  other  nonlinear  problems. 


2.  Preliminary  results.  For  later  developments,  it  is  convenient  to  recall 
certain  results  concerning  the  geometry  of  the  shell  which  here  will  be 
referred  to  its  deformed  state.  We  will  assume  at  the  outset  that  the 
coordinate  system  is  normal  relative  to  the  mid  die  surface  of  the  shell 
and  hence  use  freely  such  simplifications  as  evolve  (see,  Synge  and  Schild 
[191(9],  Sec.  2.6  and  Naghdi  [I962],  Sec.  3). 

3 

Let  h  denote  the  thickness  of  the  shell  and  x  be  the  coordinate 
along  the  normal  to  the  middle  surface.  Further,  with  reference  to 
Euclidean  3-space,  let  -  (j^  ,  j^3)  be  the  base  vectors  of  the  shell 
space  (defined  by  its  middle  surface  and  -  ~  S  x  ^  2  ^  ia 
its  middle  surface.  Then 


a 


(2.1) 


where 


x%« 


(2.2) 


b^„  is  the  second  fundamental  tensor  of  the  surface  and  is  the  Kro- 

P  P 

necker  delta. 

In  all  that  follows,  we  assume 


<  R 


min 


(2.3)* 


For  details  of  this  restriction  on  the  space  of  normal  coordinates 
see  Naghdi  [I962]. 
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with  denoting  the  leeet  redius  of  curveture  to  the  middle  surfece.  It 


..a 


then  follows  Inmedletely  that  the  characteristic  roots  of  x  b  are  less 

“1  ^ 

than  1  and  hence  (see  Mlrsky  [19^3]  P8«  332)  a  Is  given  by  the  con- 

P 

vergent  series 


- 1  ^  ■ 


P  -  0 


(2.U) 


a 

where  the  notetlon  b  ^  stands  for  the  P  factors 

P 


b“  -  b“  , 
P  v’ 


(2.5) 


Here,  we  also  note  that  the  determinant  |i  of  |i  Is  given  by 

P 


1  -  2  H  x^  +  K(x^)^  , 


2H 


uOt 


a 


and  that  If  T^,  Is  an  arbitrary  space  tensor  with  subtensor  lP  ,  then 
J  ^ 

the  corresponding  surface  tensor  r  ^  Is  defined  through  (see  Naghdi  [1962], 

Sac.  3.3). 


(2.7) 


which  characterises  the  transformation  relation  between  space  and  surface 
tensors. 

He  recall  that  the  constitutive  equations  for  an  elastic  body  In  the 
sense  of  Green,  l.e.. 


a 


tj 


(2.8) 


for  an  Isotropic  material  may  be  shown  to  have  the  form  (see  e.g.. 
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Truesdell  «nd  Toupln  [1960]»  p.  730) 


0^^  -  ♦o  +  ^2  ,  (2.9) 

where  end  e^j  ere,  respectively,  the  stress  end  the  streln  tensors, 

£  is  the  streln  energy  and  In  (2.9)  (A  >  0,l,2)i  are  functions  of 
three  Independent  Invariants  of  e^^.  For  the  purposes  of  this  paper, 
the  most  convenient  set  of  Invariants  Is  the  set 

-  e^^  ,  Ka  -  I  e^^  e^^  ,  Kg  -  ^  e^^  e^j^  e*^^  .  (2.10) 

As  stated  previously.  It  Is  the  purpose  of  the  present  paper  to  derive 
general  constitutive  equations  wherein  such  nonllnearltles  as  exist  are 
purely  physical  and  are  cast  Into  the  coefficients  of  (2.9) •  Accordingly, 
we  confine  attention  to  linearized  straln-dlsplacenent  relations 


^  ‘IJ  ’  “1|J  '"“Jll 

and  In  keeping  with  this  plan  set  ^2  ”  ^  (2.9)  •  then  follows  from 

(A. 6)  of  the  Appendix  that  £  Is  Independent  of  Kg  so  that  the  relevent 
Invariants  are 


-  »“la  *  “'I,  - 
i  t  ““Ip 

♦“•“a  +  . 

where  u^  are  the  components  of  dlsplecements  and  (  )|^  denotes  co- 

variant  differentiation. 

The  analysis  that  follows  Is  developed  under  the  Love-Klrchhoff 
hypothesis,  although  It  suy  also  be  csrrled  out  so  as  to  Incorporate  the 
effect  of  transverse  shear.  It  has  been  shown  (Naghdl  [I962])  that 
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under  the  Love -Kir chhoff  hypotheele  the  uiuel  dleplecement  aesumptlone  In 
the  form 


u  -  v(x,x)+x3(x,x). 


3  .12. 

U  -  W  (x  ,  X  ) 


(2.13) 


are  in  fact  exact.  With  (2.I3)  apd  the  use  of  (2.11),  we  have 

(  r%  *  .%v  • 


2  ‘oe  ‘  “‘’o  (  'oB  *  **  '^ps’  *  “"b  <  V  *  V  >  • 


(2.14) 


and 


Sp  ■  *03  -  0  . 


Ot  Ctii  .oe 

1’ B  • "b  •  “  "  B  • 


ft 


a 


hla  +  ^a  V  • 


(2.15) 


where  (  )||^  denotes  covarlant  derivative  with  respect  to  the  first 

fundamental  tensor  of  the  surface  a 

op 

In  view  of  (2.13),  the  invariants  and  Kg  take  the  form 


ia 


_  -la  -16  r  3  .^P  3  i  ®p  , 

Ka  -  ue  p  (  ^ap  +  *  ®ap  ^  “op  ^  * 


(2.16) 


where 
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op 

ep 


1  r  0  P  ,  ♦  ®P  1 

T;  f  p  a  +  a  ’'♦p  •  ^  • 


B. 


ap 

PP 

®ap 


i  I  /p  "‘’a  +  '<%  A  *  '.9  ♦  "*0  Ue'>  •*'’  >  •  <^-‘'^> 


1  r  ^  P  ^  Pp  1 

4  p  a  a  ♦p 


exhibit  the  symnetrles 


«p  pe 

ap  "  pet  ’ 


-«P  .  ,P® 

®OP  ®PQ! 


pp  pp 

E  «  E 

ap  pa 


(2.18) 


The  etrees-etreln  reletions.  He  have  already  admitted  the  linearized 
strain-displacement  relation  (2.11),  the  displacements  in  the  form  (2.I3) 
which  are  exact  under  the  Love-Klrchho££  hypothesis,  and  moreover  have 
required  0.  On  this  basis,  we  now  proceed  to  derive  exact  constitutive 
equations  for  elastic  shells  which  are  nonlinear  In  the  sense  stated  In 
Section  1. 

First,  we  recall  the  definition  of  the  stress  resultants  In  the 

form 


h 


where  (2.8)  has  been  used.  Observing  from  (2.I3)  that 


^ 


ap 


d  e 

pv  'op 


■  ^  5  (  “  6®  +1.“  6®  ) 


8  L  a 

•  II 


de. 


a  V 


(3-2) 
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then  n*y  be  expressed  as 


■I 


dx 


3 


L  dx 


and  If  we  define 


then 


(3.5) 


(3.»^) 


(3.5) 


The  expression  {^.3)^  Involving  the  stress  couples  follows  from  the 
saaw  procedure  used  to  obtain  (3«3)x  • 

We  now  proceed  to  obtain  a  suitable  representation  for  the  scalar 
potential  X  In  terms  of  the  Invariants  and  Kg.  Because  of  the  re¬ 

striction  ^2  •  0  ,  we  have  from  (A. 9)  that  C  ^  .  =  0.  It  then  follows  from 

Ar  A 

(A. 13)  that  L  has  the  form  ' 


Z  - 


1  _  „A  j^+1  ^  1  ,  „A+ 1 


TT  V 


A+1  "'A*'! 


Consider  now 


p  -  0 


p  -  0 


(3.6) 


(3.7) 


where  (2.4)  has  been  used.  In  the  last  series  we  replace  P  -f  1  by  P  and 
write 


P  -  1  P  -  1 
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(3.8) 


so  that  if  we  define 


P 

M 


{ 


,3 

3 

> 

for  P  -  0  , 

3 

**a 

b“ 

A  A®  V*,  . 

for  P  £  1  5 

a 

3 

a  3 

then  (3 >3)  becomes 


Ki 


P 

M 


> 


(3.9) 


(3.10) 


The  right  hand  side  of  (3.IO)  is  a  converging  series  since  the  series 

for  converges.  A  parallel  procedure  gives  for  the  expression 

3 


Kg  -  ^  (x®)^  N  , 


P  ■  0 


(3.11) 


where 


P 

N 


/  *®P 


. 

a  vp  ’ 


V 


'Qp  a3 


p  ap 


®  P  ‘03  ’ 


for  P  -  0  , 


for  P  -  1 


for  P  2  2  . 


(3.12) 


In  (3.12),  A,  B  and  E  are  given  by  (2.17)  and 

k-t  -Z  ■ 


Q  -  0 


Q  •  0 


(3.13) 


where  it  should  be  noted  that  C  is  symnetric,  i.e.. 


% 


Cpo 

'P0 


(3.H) 
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By  a  procedure  similar  to  that  used  In  the  foregoing,  it  can  be 
shown  that  the  powers  of  the  invariants  K],  and  K2  may  be  expressed  as 


■I 


/  3vP 

(x  ) 


P 

D  , 
A 


where 


S  T-S 
M  M  . 


P 

D 

A 


P  0  T 

S  T-S 

I  •  •  •  I  "  " 


Q  -  0  R 


Q  -R  P-Q 
.  M  M 


9 


Q-R  P-Q 
M  N  . 


(3.15) 


(3.16) 


The  formalism  of  (3.I6)  is  intended  to  indicate  that  there  are  A  factors 
and  A  *  1  summations.  It  also  follows  from  (3«l6)  that 


P 

C 

A 


Q  -  0 


P-Q 

M 


For  the  product  (Kj^)^  (Kj)^  ,  we  write 


(K,)^  (K^f  -  f  (xY  A  , 

_  ^  A  n 


p  -  0 


A  r 


p 

A 

A  r 


.  y  ?  'i’  .  y  8 

^  A  r  ^  r  z 


Q  -  0 


Q  -  0 


-Q 

c 

A 


(5.17) 


(3. 18) 


and  in  order  that  the  representation  (3.I8)  be  applicable  to  the  special 
cases  A  ■  0  and  r  ■  0  ,  we  append  to  (3* 16)  the  definitions 
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0  0  P  P 

C*Dal,  C>D»0,  for  P  2  1  , 

0  0  0  0 

p  p 

C-D-0,  £orP<0 
A  ^ 


(5-19) 


P 

and  note  In  passing  that  A  is  independent  of  . 

A  r 

With  the  aid  of  (3.18)  and  (2.6),  (3.4)  may  now  be  written  as 


■/ 


h  “  o 

2  ,  3  /  3v2  ,  r  1  ,  3vP  ^ 

[  1  -  2H  X  +  K(x  )  ]  [  2^  (x  y 

P  -  0 


r+1 

(3.20) 


\  .  L  y  (x^)^  A  ]  dx“ 
‘  +  1  A  A+1  0 

P  -  0 


Introducing  the  notation 

Ti  (Q)  -  i(|)‘^  (1  -  (-1)*^  ) 

and  observing  that 

Tj  (Q)  ■  0  ,  for  Q  even, 

TJ  (Q)  -  ^  (  J  ^  .  for  Q  odd  . 

then  the  integration  of  (3.20)  yields 
00 

-■If 

P  -  0  ■ 

■( 


[  T)  (P  +  1)  -  2Htj  (M-2  )  +  K  n  (P  +  5)  ] 


p  1 

^  ^  A  r  +  1  ^ 


^aJo 


(5.21) 


(3.22) 


(3.23) 
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In  order  to  express  tht  stress  and  couple  resultants  and 

0  0 

in  terms  of  the  invariants  C  and  ^  ,  we  need  to  obtain  certain  additional 
results.  To  this  end,  from  the  definitions  (3.9)  we  obtain 


0 

d  M 


0 


(3.2J+) 


P  P  -  1 

dJl -  ^ 

b  K 


(P  i  1). 


It  follows  from  and  (2.I7)  that 


0 

M 


(3.25) 


and  hence, 


0 

d  N 


^iV 


+  ). 


0 

B  N 
a  K 


-  0 


(3.26) 


It  is  convenient  to  introduce  the  quantity 


a 


b''  + 

P  P 


in  terms  of  which 


1 

N 


so  that 


1 

a  11 

1 

a  N 

a  »c 

l^v 


i  [  b''p  +  yPP)  +  )  j  ^ 


(3.27) 


(3.28) 


(3.29) 
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Similarly,  for  P  2  2,  we  have 


P 

N 


1  r  o  \  ^  V  /  0(  .  Qv  , 

T[2\,b  „(7d,+7*)  + 


+  (P  ..  1)  + 


P  -  2 

I 

Q  -  0 


\  b''  ] 

V  p  \  ® 


(3.30) 


so  that  for  P  2  2, 


P 

a  N 


S  7 


pv 


p 

a  N 


pv 


i  [  £“„  (  ♦  r>^)  +  p 

i 


2  ■  a 

P  -  1 


Q  -  0 


+  a^ 


P  -  2 

,0  Q. 


I 

Q  -  0 


P  ^ 


3^^ 


]  . 


(3.31) 


We  now  turn  to  A  in  (3.I8)  and  proceed  to  establish  the  following 

A  r 

identities: 


A  r 


A  r 


P  -  z 


a  A  /  an  -  A., A 


A-ir 


p  -  z 


a  A„  /  aw  -  r  A 


(3.32) 


A  r- 1 


To  prove  the  first  of  (3*32),  we  first  note  that 


P  P 

C  -  N  . 

1 


(3.33) 
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Hence, 


P  Z 

a  C  /  d  M  - 
1 


^'z 


P  -  Z 
C 

0 


(5.51^) 


Next,  from  (j.l?) 


P 

C 

2 


r 

■  I 


Q  P-Q 
C  M 
1 


Q  -  0 


(3.55) 


so  that 


P  Z 

a  C  /  3  M 

2 


r 

- 1 

Q  -  0 

P-Z 
-  2  M 


P-Q  Q  P-Q 
M  +  C  6^ 


P-Z 

2  C 
1 


(3.36) 


and  by  Induction  we  arrive  at 

P  Z  P-Z 

3C/dM  -  A  C 
A  A  -  1 


(3.37) 


Now,  since  by  (3.I8) 


P 

A 

A  r 


r 

F 


Q  -  0 


Q  P-Q 
C  D 
A  r 


(3.38) 


and  since  p  as  given  by  {^.l6)  is  dependent  only  on  the  N's,  we  have  in 
view  of 


Q  -  Z  P-Q 


P  Z  ^ 

3  A  /  a  M  -  ) 

^  A  -  1  r 

Q  -  0  ^  ^ 


A  C 


A  r 


(3.39) 


If  we  replace  Q  by  Q  4-  Z,  (^.39)  becomes 
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P  z 

S  /  S  M 

A  r 


A 


i 


Q  P-Z-Q 
C  D 

A  - 1  r 


and  then,  through  the  definitions  (3-l9)»  the  first  of  (3.32)  follows. 
The  second  of  (3.3^)  may  be  established  In  a  similar  manner  where  In 
particular  It  Is  convenient  to  use  (3.18)3, 

Explicit  forms  for  the  constitutive  equations  can  nov  be  obtained 
from  (3.5)  and  (3-23)  and  read 


00  P 


^  ^  Itllfw-l)  -  2Ht,(P«)  +  Kn(P+5)l  »  /  8r| 


P  -  0  Z-0 


pa 

(5.W) 


■  —  V  /[  ^  .IL  a.  V  ^  ^  ^  1  \ 

r  +  1  A-l  Vl  ^  V  a\  by^  ^  f 


00  p 


y  Y  |[  Ti(P+l)  -  2Hq(P+2)  +  Kti(P+5)]  | 


P -0  Z-0 


P-Z 


r  +  1  A-^iM-i  a^ 


8 

^pa 


p-  z 


i1 

'■  ^  J 


Z 

b  M 


pa 


(3.»(2) 


The  constitutive  relations  (3.UI)  and  (3.I42)  may  be  put  In  more  explicit 
forms  by  further  substitution  from  results  of  the  type  (3.2h),  (3.26), 
etc.,  but  we  delay  this  until  the  next  section,  where  special  cases  are 
considered. 

Before  closing  this  section,  we  examine  the  invariance  properties  of 
the  above  constitutive  equations.^  It  is  evident  that  the  equations  (3.41) 
and  (3.42)  are  tensorlally  Invariant.  Furthermore,  It  is  easily  verified 
that  they  are  also  consistent  with  the  virtual  work  theorem  for  shells. 
Next,  It  Is  necessary  to  show  that  the  sixth  equation  of  equilibrium. 


For  an  account  of  Invariance  requirements  in  shell  theory,  see 
Naghdl  [1962]. 
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namely 


ap  ' 


1,“  *'«) 
V 


0  , 


(3.»*3) 


or  equivalently 


ap 


(5. 41*) 


la  Identically  satisfied.  To  this  end,  observing  that  X  may  be  regarded 
as  a  function  of  &  and  N  ,  by  (3.14-4)  we  have 


00 

Z 

z  -  0 


ap 


-  b“ 


.  Z  Z  z  ^ 

ax/SM  )  +  (  ax/SN) 

^  °  'pot  ®  'pa 

[(3X/SH)  ^  (ax/aJ)  1^]} 


(3.  *^5) 


Z  -  0 


%  ( 

Z 

'  an 

.  kO.  |i_\ 

z 

z 

z 

+  d  X/ d  N 

S  N 

V  ^  ^pa 

■  Sv)} 

Consider  now  the  coefficient  of  dX/  d  M  in  (3-^3)  conjunction  with 
equations  (3.24).  For  Z2  1,  it  is  identically  aero  while  for  Z  ■  0,  it 
is  syomwtric  in  a  and  p  and  hence  vanishes.  Next,  consider  the  coefficient 
of  d  X  /  d  N.  A  little  manipulation  of  (3>3l)  will  show  that 


Z 
d  N 


-  b' 


2 

V  Tk  o  ^  °  a  ^ 


pa 


(3.>*6) 


pv 


an  expression  symmetric  in  a  and  p.  Thus,  for  Z  2  2,  in  view  of  the 

presence  of  e  in  this  coefficient  also  vanishes.  The  vanishing 

ap 

of  this  coefficient  for  Z  -  0.,  1,  follows  with  the  aid  of  equationa  (3.26) 
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and  (^.29).  The  final  Invariance  requirement  which  should  be  examined  is 
invariance  under  rigid  rotation.  It  may  be  shown  that  the  value  of  each 
of  the  invariants  M  and  N  is  unchanged  by  infinitesimal  rigid  rotations  and 
thus  X  retains  its  value  under  such  a  transformation.  The  calculations 
are  somewhat  lengthy  and  hence  are  not  included  here. 


4.  Approximate  constitutive  equations.  From  a  practical  point  of  view,  it 

is  of  interest  to  consider  constitutive  equations  which  are  approximations  to 
those  given  by  (3* 41)  and  (3.42).  The  procedure  to  be  used  in  obtaining  the 
approximate  equations  will  be  one  of  simply  terminating  the  series  after  a 
finite  number  of  terms  and  we  proceed  to  do  this  in  a  systematic  manner. 

First,  we  recall  that  there  are  two  types  of  summations  involved  in 
(3.41)  and  (3.42),  summation  over  the  letter  P  and  summation  (implied) 
over  the  indices  A  and  r.  These  two  summations  are  effected  independently 
and  it  is  this  fact  which  leads  to  the  independent  consideration  of  geo¬ 
metric  properties  and  material  response.  In  order  to  systematize  our  dis¬ 
cussion  of  this  matter,  we  will  use  the  notation  —  to  represent  terms  of 

CK  ”  d 

the  type  b  and  7  to  represent  terms  of  the  type  7  .  Terms  of  the  type 

QL  OL  ^  y  ^ 

K  .  and  \  will  be  regarded  as  j.  If  we  return  to  the  definitions  of 

P  P  K 

the  various  elements  which  are  the  constituents  of  (3.41)  and  (3.b2),  we 
can  easily  see  to  what  power  7  and  ^  occur  in  them.  These  results  are 
summarized  in  Table  I  below: 


Element 


P 

M 

P 

N 

P 

C 

A 

P 

D 

r 

p 

A 

A  r 


TABLE  I 

Power  to  which 
7  occurs 

1 

2 

A 

2r 

A  +  2r 


Pow^r  to  which 
—  occurs 

P 

P 

P 

P 

P 


-17- 


We  observe  from  this  table  that  the  powers  to  which  y  and  ~  occur  In 

the  various  elements  are  Independent  In  the  sense  that  the  former  Is  determined 

by  A  *nd  r  while  the  latter  Is  determined  by  P. 

The  quantities  and  In  (3*^1)  ond  (3*^)  material  coefficients, 

and  the  power  to  which  y  occurs  in  terms  involving  these  coefficients  Is 
determined  solely  by  the  value  of  the  Index  of  those  coefficients.  This 
fact  is  seen  by  applying  the  information  of  Table  I  to  (3AI)  and  (3«te) 
and  the  result  is  summarized  In  Table  II  which,  together  with  Table  I, 
reveals  the  manner  In  which  the  y  and  ^  terms  occur  In  the  constitutive 
equations. 

TABLE  II 


Coefficient  of 


Power  to  which 
y  occurs 


A 


AT 


A  +  ar  +  1 


With  the  Information  gleaned  from  the  above  two  tables,  the  pro* 
cedure  for  determining  various  approximate  constitutive  equations  Is  now 
clear.  First,  we  select  the  highest  power  to  which  we  wish  —  to  occur  in 
the  expansion  and  we  then  sum  over  P  in  such  a  manner  as  to  obtain  all 
terms  up  to  and  including  that  power.  We  emphasize,  however,  that  because 
of  the  presence  of  q  (P  -f  3)  in  (3.UI)  and  (3.I42),  we  cannot  choose  the 
highest  power  of  ^  to  be  less  than  2.  After  the  choice  on  ~  has  been 
made,  we  are  free  to  select  the  power  to  which  y  shall  appear  in  the  con- 
stltutlve  equations.  The  term  "order"  will  be  used  In  this  connection, 

1  e.,  equations  linear  In  7  will  be  called  first  order,  those  containing 
both  linear  terms  and  terms  of  second  degree  in  7  will  be  called  second 
order,  etc.  It  Is  clear  that  once  the  choices  on  —  and  7  have  been  made, 

K 

the  determination  of  the  specific  equations  Is  a  routine  matter. 


The  sequence  of  choices  used  here  In  selecting  the  powers  of  ^  and  7 
may  be  Interchanged.  The  one  adopted  here  appears  to  be  convenient. 
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In  the  examples  which  follow,  we  shall  always  choose  the  highest 
power  to  which  ^  occurs  to  be  2.  Consequently,  we  will  first  expand  the 
relations  (3.41)  and  (3-42)  Imposing  this  restriction  and  then,  by  making 
use  of  relations  of  the  type  (3*46),  obtain  the  desired  expressions.  Thus 
expansion  of  (3.41)  In  the  manner  Just  described  gives 


[r  t  1  A-i^r+1 


1 


5^  Vo  FF 


0  V  S 

pa  A  0  o  7pQ( 


1  0  0  1 

.2H(A,  A  ^ 

AO  ^7^  AO  ^7^ 


)  +  K  A 
A 


..ULl 

0  ^  V  -1 


r  +  1  Ja- 


0  1 


1 

Am 


aL 


I'WI  a  ’aVr+1  5^  *i.i*r+i 


•  k-fr+i  ^  ^  >  *  ViVi  ^  ] 


+  F, 


r  A 

^r\  A  r 


0 

a  H 


1 

+  A 


1 

AJL 


A  r  ^  7| 


pa 


0 

+  ,A 
A  r 


^  r, 


pa 


1  0  Q  1  0 

.  2H  (  A  ^  ^  ? 

6  r  51^  *  r  i  r 


The  Cayley-Hamllton 
may  be  expressed  as 


theorem  when 


applied  to  b^’  with  the  use  of  (3.24) 
P 


(4.2) 


and  If  we  write 


p  for  the  adjoint  matrix  of  b^^,  we  have 
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(It- 5) 


0  1 
.  2,^  .  ^ 

^  ’’so  ^  V 


Further,  combination  of  (4.3)>  (3*29)  >nd  (3>3l-)  yields 


2H 


1 

N 


+  K 


0 

N 


pot 


.  I  [  bP”  -b  b®^  b“''  - (b''»  -f  1 


(ii.it) 


end 


-i  (x.“e  5?^  (rP''b/»)l  (b.5) 

.  i  ((b“^  /'■  *  b9^  r)  -  (r®''  +  r^)  *  (<“»  +«»)). 


Now,  with  the  use  of  (4.2)  to  (4-3)  end  with  some  simplification  and 
combination;  (4.1)  reduces  to 


L*o  •“*'  a'V)  ®“^]  *  rVi  'a^[a-^+l  A-i*r-n 

p«] 

■  1,  (r“®  +  r®“)  ♦  -  tP,(r'«  +  r®'') 

(It.  6) 

) 

( 

'(b“^  b  bS^  b“'')  -  (b'lS  +  x?'’))  ]  }  . 

-20- 


and,  similarly,  the  corresponding  expression  for  is 

^FtT'.rUv. 

*  f  v[a  ^  -  «“v  *  y'^'^)] }  • 

Equations  (4.6)  and  (4.7)  are  of  second  degree  in  ^  but  of  arbitrary 
order  in  y.  Such  subsequent  examples  as  are  to  be  given  will  follow  directly 
from  these  equations  after  a  choice  has  been  made  on  the  order  of  7.  As 
a  first  example,  let  us  obtain  the  first  order  or  linear  constitutive  equa¬ 
tions.  For  this  case.  Table  II  instructs  us  to  retain  Lj.  and  and  to 
set  all  other  L^,  p  equal  to  aero.  (In  order  that  the  state  of  zero 
stress  and  zero  strain  coincide,  we  also  set  ■  0.)  In  fact,  if  we  set 


Ll 


1  -  V® 


f 


'00 


E 

1  +  V  ’ 


(4.8) 


we  then  obtain  the  Flugge-Lur  'e-Byme  constitutive  equations.  For  future 
reference,  we  record  representative  examples  of  these  equations,  written 
in  lines  of  curvature  coordinates  and  in  physical  components.  Thus, 

•*(11)  ■  ’'(11)  ^1^(22)  12  Ri^^'^(ii)  "  Ri’'(li)^  ’ 

”(12)  "  2  ^00  {^’'(i2)'’’’'(ai)^'''  12  ^Rg  ’  Ri^^'^(2i)  ■  Ri  ’'(21]^} 

(4.9) 

•*(11)  “  I2  {^^^■'■^00^  '^(11)  ^i'^(22)  ^^i'*’^00^^Rg  ■  Ri^  ’'(11]^ 

“(12)  ■  2)1  ^00  {  ^*^(12)  '^(ai))'*’  (r2  ■  Rj^  ^(21)  } 

where  the  constants  and  F^^  are  given  by  (4.8). 

We  now  return  to  equations  (4.6)  and  (4.7)  and  proceed  to  obtain  the 
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constitutive  equations  for  the  second  order  theory.  For  this  purpose, 
it  is  convenient  to  introduce  the  notation 


(4.10) 


where  the  subscript  k  in  each  case  designates  that  resultant  which  arises 
from  (4.6)  and  (4.7)  ^rom  exactly  those  terms  which  are  of  degree  k  in  y. 
Thus,  and  would  refer  to  the  right  hand  sides  of  the  Pltlgge- 

Lur'e~Byme  equations  (4.9)*  If  t4e  theory  sought  is  order,  then  we 
have 


k  -  1  k  -  1 


while  for  a  completely  general  theory. 


OB 

■I 

k  -  1 


(4.11) 


(4.12) 


With  this  understanding,  we  need  only  determine  and  to  complete 
the  description  of  the  second  order  theory  since  and  x*P^  given 
4y  (4.9).  The  representative  results  are  as  follows: 

12  {*'*^^^’'(11)  ’'(2^  ^  ^(Xl)  ■  Ri^  ^^(ix)  ^22)^  ^ 

^10  [^Rg  "  Rx^**  ^(IX)  ’'(Xl)  ^  ^(21)  ^’'(X2)  ’'(21)^  )  (4.13) 

*  ^Ra  ■  Rx^  Nxx)  ^^’'(xx)  +  ’'(22)^  +  ^^(XX)  +  ^(22)^  ^(n)]  }  - 
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2**(i2)  "  2  ^10  '{  ^(22)^  ■*■  ^(21)^ 

(»V.U) 

+  12  [^\ll)  ■^^22)^^^12)  +^2l)^  +  ^’'(11)  +  ’'(22)^  ^2i)  ‘  }  ’ 

and 

2«(11)  “  fe{^2(7(ii)+  7(23))  [^^Nll)  ^22))+  i  ^’'(n)  +  ’'(22))] 

+  ^0  [^\  !!)’'( 11)  ■^\22)’'(22))+  2  ^\l2)  +^21))  ^’'(12)  '^’'(21)) 

(1^.15) 

+  ■^’'^2))  +^(’'(12)  +7(21))®) 

+  ^\ll)  +  ^22))  ’'(H)  +  <’'(11)  +  ’'(22))  ('^(11)  +  ’'(ii)^R2  ■  i^)]  }  ‘ 

^(X2)  “  24^10^^11)  *^^22))  ^’'(X2)  ■^’'(2X)) 

(U.16) 

+  ^’'(Xl)  +  ’'(22))  ['"(is)  +  '"(21)  +  ’'(2X)  i  )  ]  }  • 

The  foregoing  expressions  simplify  somewhat  for  special  configurations.  For 
example,  in  the  case  of  a  flat  plate,  the  representative  equations  reduce  to 

2*(XX)  “  **  {*'2^’'(tl)‘’‘  ’'(22))  >^10  [2  ^’'(Xl)  ’'^22)) 

+  I  <’'(X2)  +  ’'(2X))^  +  ’'(H)  ^’'(XX)  ’'(22))]  } 

(1^.17) 

+  ^|l2  +  *^(22))^  *  ^10  [^^'^(11)  '^izz9 

+  (*^(12)  +'^(2X))*+''(1X)  ('^(XX)  +'^(22))]  }  » 
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(iz) 


2  ^10 


^(22)^  ^^(12)  "*■  ^(sx)^ 


(4.18) 

+  ^<*(11)  +‘(2,))  <'(12)  +'‘(21))}  ■ 

^(11)  •  ’'(22)*<'‘(ll)+  “(22))  +  'ic[5  *(11)  )’(11)+  ''(22)’'(22) 

(4.19) 

+  (*(12)  +*  (21))  (>'(12)  ♦  >'(21))  +  >'(11)  *(22)  ■*■  ^(22)  *(11)]  }  • 


r 

^(12)  "  ^  ^10  *^(22)^  ^^(12)  ^  ^(si)^ 

♦('(11,  ♦>'(.»))  (*(12)  **(21))}  • 


(4.20) 


Our  purpose  in  recording  (4.17)  to  (4.20)  is  to  observe  the  coupling,  l.e., 
Che  dependence  of  the  stress  resultants  on  the  strsins  K  as  well  as  the 

Ofp 

dependence  of  the  stress  couples  on  the  strains  In  contrast,  it  aey 

be  recalled  that  when  (4.9)  !•  specialised  to  the  case  of  a  flat  plata,  no 
such  coupling  is  present. 

It  is  of  interest  to  compare  the  constitutive  equations  (j.4l)  and 

(^.42)  upon  specialization  with  those  obtained  by  Zema  [i960]  whose  con- 

## 

stitutive  equations  contain  only  first  and  third  degree  terms.  If  in 
the  expressions  for  and  iP^,  we  neglect,  as  Zerna  does,  terms  involving 
vlbh  respect  to  the  resiaining  ones  as  well  as  terms  of  the  type 
then  (4.6)  and  (4.7)  becoaie 


■  h 


+ 


AF 


(4.21) 


Although  Zema  [i960]  dotis  not  explicitly  invoke  the  Love-Kirchhoff 
hypothesis,  he  neglects  the  effect  of  7  and  sets  0^^  m  0  but  retains 
the  effect  of  733. 
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and 


^  {K  +  'Ar  4  -iV  +1  ] 

}  . 


+  'Ar  A*r  ^  ♦  'af  A°r 


ik.22) 


where  we  have  written 


(ll.25) 


Now  If,  as  In  Zerna's  case,  we  wish  to  retain  only  the  first  and  third  degree 
terms.  It  follows  from  Table  II  that  the  only  nonvanishing  coefficients  are 
1^  ,  L3  ,  Fqq  ,  FgQ  ,  and  F^^.  Here  expansion  of  (4.21)  gives 

*  [>-3(5^2)’  ♦  'a,  A>]  •“" 

'A)]  ^ 

which  Is  of  the  same  form  as  the  corresponding  expression  given  by  Zerna. 

It  should  be  noted,  however,  that  while  the  coefficients  of  {7^q)(,7^  7^p) 
and  (7  q)^  7^^  Identical  In  (4.24),  they  are  not  the  same  In  the 
corresponding  expression  given  by  Zerna  [i960].  A  similar  conclusion  can 
be  arrived  at  upon  the  comparison  of  (4.22),  when  only  first  and  third  degree 
terms  are  retained,  with  Zerna's  expression  for  the  stress  couples  after  a 
misprint  In  Zerna's  equation  (4^)  has  been  corrected. 

5.  Membrane  theory .  In  this  section,  we  specialize  the  previous  results  to 

those  appropriate  for  the  membrance  theory  of  shells.  We  recall  that 
for  equilibrium  problems  of  the  membrane  theory,  where  the  shear  stress 
resultants  and  the  stress  couples  are  absent,  the  state  of  stress  Is  stati¬ 
cally  determinate  and  the  stress  resultants  are  obtained  from  the 
equations 
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(5-1) 


lU  +  P®  ■  0  . 

'■0.9''^®  ■  °  ■ 


For  the  subsequent  determlnetlon  of  the  dlsplscements,  we  turn  to 
the  constitutive  equstlons  (l4'.6)  and,  consistent  with  the  membrsne  theory, 
neglect  the  terms  with  the  factor  (-)  Thus,  for  the  membrane  theory, 
instead  of  (4.6),  we  have 


r  +  1  A‘i^r+1^  “ 


AP  A 


(5.2) 


Because  of  the  symmetry  of  the  right-hand  side  of  it  is  seen 

that  (^-Os  identically  satisfied.  Also,  in  view  01^3.18)2,  all 
A's  appearing  in  (3*2)  may  be  written  in  terms  of  M  and  H  alone  so  that 

iPP  -  h  F^j,  (N)^  ]  a“P 

0  0  ^5.5) 

+  (N)^  +  • 

0  0 

Since  N  and  N  are  functions  of  7  it  is  clear  that  the  combination 
(3*3)  (2.11)  constitutes  a  system  of  nonlinear  differential 

equations  in  the  displacements.  Since  iP^  as  solutions  of  (3.I)  can  be 

0  0  rtfl 

regarded  as  known,  it  is  convenient  to  express  M  and  M  in  terms  of  li^, 
or  sK>re  specifically,  in  terms  of  the  invariants  of  iP^.  Let 

^  ’  ^2  •  2  **^0  *  ^5.4) 


and  introduce  further  the  notation 


AP 


(5-5) 
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Contraction  on  (5«3)  gives 

^  ■’T+T  ^Ar  ^ 


(5.6) 


or 


i  (J^-»H) 


(5.7) 


in  terms  of  which  becomes 


rf'e  -  h(i(  Si-  *2)  .“S  *  ,  i  (,“e  ^  ,eo.,  I 


(5.8) 


By  (5.4)  and  (3.8),  we  arrive  at  the  additional  relation 

4  P2  -  -  h^  /  [UH-(M)®].  (5.9) 

0  0 

Equations  (3.6)  and  (3*9)  give  M  and  N  in  terms  of  Pi  and  P.  .  In 

0 

the  linear  case,^the  solution  for  M  is  straightforward,  and  since  in  this 
case  ■  Fqq  ,  M  is  determined  by 


h  (l^L^  +  Fqq)  M  -  Pi  - 


(5.10) 


It  then  follows  from  (3.3)>  after  the  values  of  F^  and  are  supplied 
from  (4.8),  that 


5  (r“^  *7^)  -  s  ‘<1  +  V)  -  V  y,  .“P  J 


(5.11) 


which  is  the  usual  result  of  the  linear  theory. 
In  the  second  order  theory. 


^  *  F  +  P  M 

^  00  10 


(5.12) 


so  that  by  (3.6)  and  (3.9)>  «nd  Pg  reduce  to 
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(5.13) 


ro  0  0  001 

Pi  -  h  .(2  [LiM  +  L2(M)  +  N]  +  (Fqq  +  P^qM)  M  J- 


P2  -  h  (Pqq  +  PjQ  M)^  [  IrtJ  -  (M)®]  . 


(5.U) 


0  0 
Although  N  may  be  easily  eliminated  between  these  equations,  giving  M 


in  terms  of  P^^and  P2,  the  resulting  equation  ^.s  of  the  V 


th 


degree  in  N. 

However,  once  M  is  known  in  terms  of  Pi  and  Pg,  (3.3)  yields  a  system  of 
linear  differential  equations  for  the  displacements.  In  fact,  ^ven  ig 
the  general  case,  so  long  as  (3.6)  and  (3>9)  on  be  solved  for M  and  N  in 
terms  of  Pi  and  P2  ,  the  differential  equations  for  the  displacements  are 
linear. 


6.  Appendix:  Power  aeries  forms  for  the  elastic  potential. 

Equations  (2.9)  include  the  most  general  isotropic  stress-strain 
relations  of  elasticity,  but  since  the  being  arbitrary  functions  of 
the  invariants  of  e^j  depend  also  on  x®,  it  is  desirable  to  express  £  in 
(2.8)  in  a  form  which  will  readily  lend  itself  to  integration  across  the 
thickness  of  the  shell.  The  developments  here,  however,  are  entirely 
general  in  that  they  are  applicable  whenever  two  sysnetric  second  order 
tensors  are  related  through  a  potential. 

Let  the  syiametric  tensors  under  discussion  be  Q^j  end  A^^,  where  the 
first  is  an  analytic  isotropic  function  of  the  second.  He  further  essuaw 
that  a  potential  r  (A  j)  exists  such  that 


(A.l) 


Let 


Ki 


(A. 2) 


♦  '  ■ 

The  necessary  and  sufficient  conditions  for  the  existence  of  a  potential  r 
may  be  found  in  Truesdell's  [1952]  paper.  Truesdell's  conditions,  however, 
are  given  in  terms  of  a  set  of  invariants  which  are  different  from  those 
utilized  in  this  paper. 
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be  selected  as  the  three  Independent  invariants,  then 


r  -  r  (Ki,  Ks.  K3  ). 


(A. 5) 


Our  choice  of  the  Invariants  (A. 2)  Is  motivated  by  the  fact  that  their 
derivatives  have  the  form 


3. 

1 


a 

TIT, 


(A.  4) 


and  hence  (A.l)  may  be  written  as 


C,  8  ,  +  Co  A 


+  C3 


(A. 5) 


where  the  coefficients 


a  r 


(A  -  1,  2,  3) 


(A.6) 


are  power  series  In  the  invariants  K^.  Since  the  existence  of  a  poten¬ 
tial  r  has  been  aasumed,  It  follows  that 


or  equivalently 

^  Ca  _  a  Cg 

aKg 


(A.7) 


(A.8) 


which  are  In  fact  the  necessary  and  sufficient  conditions  for  (A.l)  In 
terms  of  K^. 

Since  the  may  be  written  as  power  series  In  terms  of  K^,  we  write 
C3  In  the  form 

C3  -  Ki^  k/  K3^  ,  (A. 9) 
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where  At  F,  A  ■  1,  2,  39  •  ■  •  sunmetlon  is  Intended  even  In  those  terms 

(which  will  erlse  leter)  In  which  an  index  appears  more  then  twice.  With 
Cs  given  by  (A.9)i  ws  proceed  to  determine  Ci  end  with  the  use  of  (A. 8). 
Thus 


d  C2 
8  K3 


^  C3 

a  Ka 


(A. 10) 


so  that 


Ca 


-T- 
A  +  1 


K 


r-1 

2 


A+1 


(A. 11) 


where  the  second  term  on  the  right-hand  side  represents  the  arbitrary 
function  of  integration.  Similarly,  with  the  two  remaining  conditions  of 
(a. 8)  and  with  (A. 9)  and  (A. 11),  we  obtain 


-  rir  <krA  •'9'''"“ 


.  A  -  ,,A"1  ,,r+l  ..  yA 

r  +  1  'af  *^1  *^2  hi 


Having  determined  ,  we  compute  F  from  (A. 6)  and  obtain 


(A. 12) 


F  - 


A  +  1  'AF  A 


^  A  „  F  _  A+1 
Xi  Xa  ^3 


F  +  1  AF 


-  „  A  „  F+1  .  1  ,  « A+1 

^.r  Ka  ^  A  +  1 


(A. 13) 


where  the  arbitrary  constant  of  integration  has  been  set  equal  to  aero 
without  loss  in  generality.  The  choice  of  the  set  of  invariants  (A. 2) 
and  the  starting  point  (A. 9)  was  ontivated  by  the  desire  to  depress  the 
coefficient  Cs.  However,  as  should  be  apparent,  the  procedure  is  the 
same  for  any  choice  of  independent  invariants  such  as 


I 


II 


1 

2 


III 


imp 
J  nk 


I  • 


(A.  14) 
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If,  Instead  of  (A~^),  (A. 14)  is  taken  as  Independent,  then  different 
(although  equivalent)  integrability  conditions  will  result  (Truesdell 
[19^2,  p.  I33I).  In  this  case,  the  potential  has  the  same  form  as  (A.I3) 
except  with  different  coefficients  while  the  co¬ 

efficients  of  (a. 3)  are  [compare  with  (A. 11)  and  (A. 12)] 


r  -  r  TT^  ^  TT1 

A  7  TA-l„r_„A  +  l 

+ 

L  T* 

r+T  ‘^afa  ^ 

A  +  1  ^AFA 

III  +  F  I^'*’^  II^ 

AF 

+  ■ '  F 

^  r  +  1  AF 

^-1 

+L^  I^, 

(A. 15) 

■  -  'arA 

11/ 

■*  A  +  1  *kFA  ^ 

-  F  1^ 

*Ar  ^  ’ 

'3  ■  'arA  • 
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